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We present a preliminary result for second moment of the light cone wave function of the pion. This parameter 
is the subject of a discrepancy between theoretical predictions (coming from lattice and sum rules) and a recent 
experimental result (that remarkably agrees with purely perturbative predictions). In this work we exploit lattice 
hypercubic symmetries to remove power divergences and, moreover, implement a full 1-loop matching for all the 
contributing operators. 


1. INTRODUCTION 

The light cone wave function of the pion, 
(f>(x,Q 2 ) is the probability amplitude of finding 
a parton of a pion (moving with momentum p in 
a light cone frame) with parallel momentum equal 
to xp and transverse momentum less than Q. 
This wave function incorporates non-perturbative 
physics and plays an important role in exclusive 
hard scattering processes and in non-leptonic de¬ 
cays of heavy mesons. One example of applica¬ 
tion is the electromagnetic form factor of the pion, 
F(q 2 ), defined by 

( 7 r ( p ')| r ) V 9 K ( p )> = f((j>-p') 2 )(p + p')h ( i ) 

This form factor can, in fact, be written in terms 
of <f> and Th (the perturbative scattering ampli¬ 
tude for the constituents) as 

F(Q 2 ) =/ 4> ] {x, Q 2 )T h (x, y, Q 2 )(j){y, Q 2 )dxdy (2) 

More formally the light cone wave function can 
be defined as 01 

= F.T.(O\T{C( Zl ),q b 0 (z 2 )}\n) (3) 

where F.T. indicates a Fourier transform on z\ 
and Z 2 assuming 
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• the sum of the parallel components of the 
momenta of the two partons is equal to the 
total pion momentum p. 

• the transverse components have been inte¬ 
grated out up to momentum Q. 

The n-th moment of <f) is defined as 

<n= [ OKLQ 2 )d£ (4) 

J o 

In a typical lattice determination of such quan¬ 
tity the cut-off is provided by the lattice spacing, 
Q ~ a -1 . 

We finally wish to remark that (£°) = 1 is 
fixed by a normalization condition and (£ ) = 0 
because the wavefunction is symmetric under G- 
parity. Therefore (£ 2 ) is the first non trivial mo¬ 
ment of (j). 

This second moment can be related to 
(0|CVp|7r(p)) = fn(£ 2 )p»p v Pp + (5) 

where the ellipsis indicates divergent terms, 

O^n/p — QVpT5 DvDp q (6) 

and 

qD q' d = qDq' + qDq' 


( 7 ) 
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2. COMPUTATION 


Our computation was carried on 154 quenched 
gauge configurations generated by the UKQCD 
collaboration using the Wilson gauge action. We 
use the Clover 0(a 2 ) improved action for the light 
quarks. 

The parameters of our computation are: 

• Volume equal to 24 3 x 48. 

• f3 = 6.2 which corresponds to an inverse 
lattice spacing of a -1 = 2.67 ± O.lOGeV. 

• k values 0.13460, 0.13510, 0.13530 (corre¬ 
sponding to pseudoscalar masses of 748, 574 
and 490MeV respectively) 

• K cr u = 0.13582 and csw= 1-61 

For technical reasons we choose not to improve 
the operators and do not smear the light quarks. 
The latter choice is motivated by the fact that the 
local axial current seems to have better superpo¬ 
sition with the pion than other smeared operators 
we tried. 

In order to calculate any moment of the pion 
light-cone wavefunction, we are required to study 
the matrix elements of lowest twist local opera¬ 
tors between pion and vacuum. In a continuum 
world these operators are classified by their rep¬ 
resentation under the group of Lorentz, parity 
and charge conjugation. In the lattice-discretized 
world the Lorentz group is broken to Hi £ 0(4), 
the hypercubic group. Hence lowest twist local 
operators can mix with higher dimensional oper¬ 
ators and introduce power divergences. 

We choose a particular combination of the lat¬ 
tice operators that: 

• transforms under an irreducible representa¬ 
tion of the hypercubic group 

• does not mix with higher dimensional oper¬ 
ators. 


In particular we choose with p ^ v ^ p ^ 

p. which transforms as (1,1) ® 8 + under Hi. 
The 8 + irreducible representation would naively 
generate a term proportional to 


Pm 


Pu+Pp 


- (e^ pa Pa? 


( 8 ) 
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Figure 1. The Plot shows R as function of t, the 
timeslice, and k. The best fits are included (for 
each value of k). 


This contribution vanishes for our matrix ele¬ 
ments because of its parity. We introduce the 
following definition: 

R = C 2°(^P) (9) 

Pip 2 C^(t,p) p=(lli0) 

where C® ( t , p) is the spatial Fourier transform 
(at momentum p) of {Q^q'y^q' (0)), O is a 
short hand notation for and A = q^q' 

is the usual axial current. For large t = xq —■> oo, 
C‘ 2 (t) asymptotically approaches 


C 2 ~ 2E^)' <°IG(0)k(p)>e" B(p)t 

Using eq. (|), eq. fj) reduces to 


7 A —_ 

^ i£2\lattice _ __^2^MS 

Zo 


( 10 ) 


(ii) 


the desired second moment (without divergences) 
times a corrective matching factor to be deter¬ 
mined perturbatively. 


3. MATCHING 

A big part of this work was the determination 
of the matching factor Z° for O = O m „ p . The 
computation was performed assuming a gluon 
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mass as regulator and on-sliell external quarks 
with non-zero momentum. All 41 relevant dia¬ 
grams were expanded in Taylor up to second or¬ 
der in p, compatibly with eq. (|9j|). 

As result of our computation we find that 

|^ = 1-518 ( 12 ) 

and the coefficients are evaluated at a s (q*) for 
q* = 2/a. We tried varying q* by factors of two 
and this gives us an estimate of the error in the 
matching of the order of 10% that we include in 
our final result. 

4. RESULTS 


The result of our calculation is 


/ ^2 \ lattice 
\S / Q—2.67GeV 

= 0.185 ±0.032 

(13) 

/f2\MS 

\S / Q—2.67GeV 

= 0.280 ± 0.049tg;^ 

(14) 


(the first error is statistical and the second is sys¬ 
tematic due to matching, quenching error is not 
included). These numbers should be compared 


with results from sum rules 

(OfW = 0.40 ±0.05 (15) 

with preceding independent lattice results 

0 ?) l Q~i C GeV = 1.37 ±0.20 | (16) 

(£ 2 >Q~lGeV = 0.25 ±0.10 | (17) 

(i 2 ) l Q%AGeV = 0.10 ±0.12 | (18) 

and the exact asymptotic value 
<£ 2 )q=oo = 0.2 (19) 


(confirmed by the Fermilab experiment E791, 
performed at Q ~ 3 — 4GeV). 

We find that our value for the second moment 
is smaller than sum rule predictions and is closer 
to the asymptotic value. 

On the one side we conclude that the present 
computation has better control of statistical and 
perturbative errors than previous computations. 
On the other side we strongly feel that dynami¬ 
cal quarks may be playing an important role in 
this process and quenching errors are yet to be 
estimated and removed. We also believe that the 


same analysis should be repeated for different val¬ 
ues of the lattice spacing in order to study the Q 2 
dependence of (£ 2 ). 
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